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Uniform shear flow is a paradigmatic example of a nonequilibrium fluid state exhibiting non- 
Newtonian behavior. It is characterized by uniform density and temperature and a linear velocity 
profile U x (y) = ay, where a is the constant shear rate. In the case of a rarefied gas, all the relevant 
physical information is represented by the one-particle velocity distribution function /(r, v) = /(V), 
with V = v — U(r), which satisfies the standard nonlinear integro-differential Boltzmann equation. 
We have studied this state for a two-dimensional gas of Maxwell molecules with a collision rate 
K{9) oc lim e ^o £~ 2 6(8 — e), where 9 is the scattering angle, in which case the nonlinear Boltzmann 
collision operator reduces to a Fokker-Planck operator. We have found analytically that for shear 
rates larger than a certain threshold value os t h — 0.3520^ (where v is an average collision frequency 
and a t \i/v is the real root of the cubic equation 64a; 3 + 16s 2 + 12x — 9 = 0) the velocity distribution 
function exhibits an algebraic high-velocity tail of the form /(V; o) ~ |V[ _4_CT ^ a '<&(<£; a), where 
ip = tanVy/V x and the angular distribution function $(</?; a) is the solution of a modified Mathieu 
equation. The enforcement of the periodicity condition <3?(v5; a) = <5>(y> + tt; a) allows one to obtain 
the exponent cr(a) as a function of the shear rate. It diverges when a — > a t h and tends to a minimum 
value <7 m in — 1.252 in the limit a — * oo. As a consequence of this power-law decay for a > a t h, 
all the velocity moments of a degree equal to or larger than 2 + a(a) are divergent. In the high- 
velocity domain the velocity distribution is highly anisotropic, with the angular distribution sharply 
concentrated around a preferred orientation angle <p(a), which rotates from if = — 7r/4, 3n/4 when 
a — > ath to tp = 0, 7r in the limit a — > oo. 
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I. INTRODUCTION 

In the investigation of the physical properties of fluids far from equilibrium, one usually focuses on the nonlinear 
dependence of the momentum and heat fluxes on the gradients of the hydrodynamic fields. The associated transport 
properties are related to the population of molecules with energies of the order of or less than the mean kinetic 
energy, so that molecules moving with velocities much larger than the thermal velocity hardly contribute to those 
properties. However, the knowledge of the high-energy population in nonequilibrium states is important not only 
from a theoretical point of view but also because that population may play a crucial role in processes such as chemical 
reactions with a high activation energy or in the controlled thermonuclear fusion of a confined hydrogen plasma. 

Of course, a general description of the high-energy population for arbitrary nonequilibrium states is not possible. 
Therefore, it is worthwhile gaining some insight by considering particular states. It can be fairly said that one of the 
most extensively studied nonequilibrium states is the so-called uniform shear flow (USF). At a macroscopic level, it 
is characterized by a constant density n, a uniform temperature T, and a linear profile of the x component of the 
flow velocity along the y direction, i.e., U x (y) — ay, a being the constant shear rate. This shear rate represents the. 
only control parameter needed to measure the departure of the USF state from equilibrium. At a microscopic level,u 
the USF is described by a solution of the Liouville equation with Lees-Edwards boundary conditionsja which can 
be seen as periodic boundary conditions in the local rest frame. These conditions assure the consistency of uniform 
shear, density, and temperature, even far from equilibrium. This state has been widely used to study rheological 
properties, such as shear thinning and viscometric effects. It must be borne in mind that, except in the linear regime, 
the USF is not equivalent to the planar Couette flow. In the latter, the shearing is produced by "realistic" walls in 
relative motion, so that the boundary conditions correspond to particles interacting with the walls rather than to 
the generalized periodic Lees-Edwards boundary conditions. In contrast to USF, boundary effects are present in the 
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Couette flow and, in addition, the shear rate, density, and temperature ace local quantities^ Far from equilibrium, 
the rheological properties of the Couette flow differ from those of the USF.Q 

In general, no rigorous theory based on first principles exists for the USF state. On the other hand, if one 
restricts oneself to the case of dilute gases, the most relevant physical information is contained in the one-particle 
velocity distribution function /(r, v, t) and then the Liouville equation or, equivalently, the BB£LKY hierarchy can be 
successfully contracted to a closed equation for /, namely the nonlinear Boltzmann equationfla A number of exact 
results have been derived from the Boltzmann equation specialized to the case of Maxwell molecules under USF. 
Almost half a century ago, Ikenberry and TruesdellOtl showed that the infinite hierarchy of moment equations could 
be recursively solved. In particular, they obtained the temporal evolution of the second-degree moments, which am 
the quantities related to the rheological properties, for arbitrary values of the shear rate. Truesdell and MuncasterEI 
analyzed the temporal evolution of the third-degree moments (for instance, the heat flux), which are expected to 
vanish for long times because of symmetry. They observed that for shear rates larger than a certain value some of 
those moments jiacreased with time, a feature they referred to as an instability in the heat flux solution. However, 
further analysis^ has proved that such an increase is actually a consequence of the viscous heating and that the 
heat flux does indeed vanish for long times when it is properly scaled with respect to the thermal velocity, so the 
apparent instability is then removed. More recently, explicit expressions for the fourth-degree moments have also 
been derived. E2fEj While the (scaled) second-degree moments remain finite for arbitrary shear rates, there exists a 
critical value a c of the shear rate, beyond which the fourth-degree moments diverge. The_analysis of this singu^ 
behavior has been extended to moments up to degree 36 and 240 for three-dimcnsionalllj and two-dimensiona. 
systems, respectively. These exact results show that the moments of an even degree k > 4 are divergent if the shear 
rate is larger than a certain fc-dependent critical value that decreases as the degree k increases. This behavior of 
the moments indicates that the distribution function preseiits-an algebraic high-velocity tail. This expectation has 
been strongly supported by direct Monte Carlo simulations JIM-J which also show that the high-velocity distribution 
function presents a strong anisotropy. To the best of our knowledge these high-velocity properties have not been so 
far confirmed at a theoretical level. The aim of this paper is to fill this gap. 

Since, due to the complexity of the collision term, the nonlinear Boltzmann equation is extremely difficult to solve, 
we consider here a simplified model of two-dimensional Maxwell molecules. In this model the collision rate is assumed 
to vanish for all the scattering angles except for a small angle corresponding to grazing collisions. This allows us to 
replace the nonlinear collision operator by a Fokker-Planck operatom that, on the other hand, preserves most of the 
general features of the original operator. By assuming a high- velocity tail of the form /(V) ~ V~ 4 ~ rT $(ip), where 
V = v — U is the peculiar velocity and ip = tan -1 V y /V x is the polar angle in velocity space, the Boltzmann equation 
yields a linear second-order ordinary differential equation for The periodicity condition on determines a as 

a function of the shear rate, thus confirming the assumed algebraic tail. The results show that the exponent <r(a) is 
a monotonically decreasing function that diverges when a approaches a threshold value ath and tends to a minimum 
value <7 m in ~ 1.252 in the limit of large shear rates. The first property means that for shear rates smaller than ath the 
distribution function decays for large velocities more rapidly than any power law (e.g., as a stretched exponential), 
while the second property implies that, even for large shear rates, all the moments of the form (V k ) with k < 2 + c m i n 
are finite. The orientation distribution in the high-velocity domain is concentrated around a preferred angle 

<p{a) that rotates counter-clockwise as the shear rate increases. Moreover, this distribution is infinitely sharp for the 
extreme values of the shear rate: lim a ^ ath $(y>) = \\S{<p — 37r/4) + S((p — 77r/4)], lim _ >00 = \[5{(p) + 6((f — ir)]. 
To our knowledge, this is the first time that such a detailed picture of the high-velocity behavior of the solution to 
the Boltzmann equation in a far from equilibrium state has been analytically described. 

The organization of this paper is as follows. The uniform shear flow state is described in Sec. 0. Some exact 
scaling properties in the case of Maxwell molecules are used to map the time-dependent state of the system onto an 
equivalent "thermostatted" state which reaches a nonequilibrium steady state for long times. The simple scattering 



model representing grazing collisions is introduced in Sec. III. The temporal evolution of the (scaled) fourth-degree 



moments is analyzed, the results showing that they diverge for shear rates equal to or larger than a certain critical 
value a c , as expected from previous analyses for more realistic scattering laws.E3ii.il The derivation of the high- velocity 
tail of the form /(V) ~ V~ A ~ a Q{<p) is worked out in Sec. [V, where the dependence of the exponent a on the shear 
rate is obtained from the periodicity condition $(</?) = <&(</? + 7r). Finally, the results are discussed in Sec. fv| 

II. UNIFORM SHEAR FLOW 

The most relevant quantity to determine the nonequilibrium properties of a dilute gas is the one-particle velocity 
distribution function f(r, -y-ii—Its time evolution is governed by the nonlinear Boltzmann equation. In the absence 
of external forces, it reads flltJ 
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= Q[f,f], (i) 

where I(g, 6) is the differential cross section (9 being the scattering angle) and we are using standard notation to denote 
the distribution function evaluated at pre- and post-collisional velocities. Of course, Eq. ([!]) must be supplemented 
with the appropriate initial and boundary conditions. 
Let us now introduce the velocity field 

Ui(r) = aijrj, a,ij = aS lx S jy , (2) 

where a is a constant shear rate. We define the uniform shear flow (USF) state as the one that is spatially homogeneous 
when the velocities of particles are referred to a Lagrangian frame moving with the velocity field U(r), i.e., 

/(r,v,t)=/(V,i), (3) 

where V = v — U(r) is the peculiar velocity. Consequently, the Boltzmann equation (|l|) becomes 

The usual boundary conditions used to generate tbi-USF are the Lees-Edwards periodic boundary conditions ,i@ but 
the so-called "bounce-back" boundary conditionstZHla are also consistent with the USF. It is worthwhile noting that 
Eq. (|J) can be interpreted as representing a homogeneous state under the action of the nonconservative external force 
Fi = —maijVj. Note also that Eq. (||) is invariant under the transformation (V x , V y ) — * (—V x ,—V y ). 

Let us now particularize to the case of Maxwell molecules, i.e., particles interacting via a potential 4>(r) oc r~ 2 ^ _1 ), 
where d is the dimensionality of the system. In that case, the collision rate gl(g,9) = K(8) is independent of the 
relative velocity gJij'Ej From a mathematical point of view, this makes the Boltzmann collision operator more tractable 
than for other interaction models. In particular, any collisional moment of degr.ee k can be expressed as a bilinear 
combination of moments of / of degrees k' and k" , such that k' + k" = kntBlEl This allows one, in principle, to 
solve reyuiaiKely the hierarchy of moment equations arising from Eq. (J4|) , even though the expliciJ—form for / is not 
known .Lro£3 In addition, Eq. (^) for Maxwell molecules exhibits the following scaling propertyJlJ'tEl Let us introduce 
the scaled quantities 

V = e- at V, (5) 



f(V,t) = e dat f(V,t), (6) 
where a is an arbitrary constant. Then, Eq. (|4|) reduces to 

-/ - -=- {a ij V j + aViff = Q [/, /] . (7) 

This equation can be interpreted as the one corresponding to USF in presence of an external drag force —maV. Since 
the mapping of Eq. (Q) onto Eq. (Q) (and vice versa) is an exact property, we are free to choose the parameter a as 
we like. For convenience, we choose a as a function of the shear rate a by the condition that the scaled temperature 

T ( /) = < - / 'ivF 2 /(v./i (m 

reaches a constant value in the long-time limit. In the above equation, m is the mass of a particle, ks is the Boltzmann 
constant, n is the number density, and T is the unsealed temperature. With this choice of a, the term -maV plays the 
role of a thermostat force. This kind of thermostat forces is usually employed in nonequilibrium molecular dynamics 
simulations .E3 c3 Henceforth, we will adopt the point of view behind Eq. (M), i.e., all the quantities will be understood 
to be scaled quantities, and we will drop the bars for convenience. Taking second-degree moments in Eq. (0) we have 

d 

—Pij + (a ik P jk + a jk P ik ) + 2aPij = -v{P^ - pSy), (9) 

where 
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Pij = mj dVV^f (10) 

is the pressure tensor, p — nksT = 4Tr P is the hydrostatic pressure, and v is an effective collision frequency defined 
as (see Appendix [A|) 

v = n-^j J dfl K(0) sin 2 t cos 2 ~. (11) 

It is convenient to choose v~ x as the time unit and introduce the dimensionless quantities t* = vt, a* = a/v, and 
a* = ajv. Henceforth, these reduced variables will be implicitly assumed and the asterisks will be omitted. 
From Eq. (^) it is easy to get the following closed differential equation for the temperature: 

i + 2a](i + 2a + l) 2 T=^T. (12) 



dt J \dt 

So far, a is arbitrary. Now, as said before, we choose a under the condition that the temperature reaches a stationary 
value in the long-time limit. This implies that a[a) is the real root of the cubic equation 

a 2 = da(l + 2a) 2 , (13) 



i.e., 



a(a) — — sinh 2 
o 



icosh 1 (l + z) 



1 / / \V3 1 / / \ 1 / 3 1 

- (l + z+ V2z + z 2 J +-(l + z- V2z + z 2 ) -g, (14) 
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where z = ^ a 2 . The stationary values of the elements of the pressure tensor are then easily obtained from Eq. 



F l + 2a' yy l + 2a v (1 + 2a) 2 



These quantities exhibit non-Newtonian effects: normal stress differences and a nonlinear relationship between the 
shear stress and the shear rate. 

While the second-degree moments are well defined for anjj— value of the shear rate a, moments of degree four and 
higher diverge if a is larger than a certain critical value a c .E3~t3 This suggests that the stationary solution to Eq. 
(0) presents a high- velocity tail of the form /(V) ~ y- d -2-o-(a)^ w here a {a) is a positive definite exponent that is a 
decreasing function of the shear rate. 

III. A SIMPLE SCATTERING MODEL 

In the Monte Carlo simulations as well as in the theoretical analysis of higher-degree moments presented in Ref. 
[l3| the scattering law was assumed to be isotropic: K(9) = const. However, this choice for the collision rate does 
not simplify the collision operator in a significant way as to allow one to confirm theoretically the high-velocity tail. 
Since the main features can be expected to be to some extent independent of the particular scattering model K(9), it 
is convenient to consider the simplest choice lending itself to an analytic or semi- analytic treatment. Here we consider 
the following simplified scattering model: 

K{6) = K Q Yiia ( r d 5{e - e), (16) 

£-+0 

which corresponds to grazing collisions of Maxwell molecules. It is proved in Appendix ^ that in this case the collision 
operator becomes 



ov l mn oVi oVj 



/(V), (17) 



where 
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£3 = (v*£) -m^-^m 



(18) 



The right-hand side of (JT^) has the structure of a Fokker-Planck operator with a velocity-dependent diffusion 
term. The model ( |l7| ) was derived (with d = 3) in Section II. 9 of Ref. || as the contribution of grazing collisions to 
the collision operator. Since here only grazing collisions are considered, the above contribution becomes the entire 
collision operator. 

Now, we restrict ourselves to the two-dimensional case (d = 2), so that Eq. (0) becomes 
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d 



1 



Slf ~ aVy^rf - a—Vif = - —Vi 



dt 



dVi 



d 



2p8 t 



Pi j d d 



dVi 



mn 



f. 



where we have already made v = 1. Of course, by taking second-degree moments we reobtain Eq. 
now the fourth-degree moments 



M s = i J dVVt s V y s f(V), 



s = 0, . 



.4. 



(19) 

Let us consider 
(20) 



Multiplying both sides of Eq. (M) by V% s V y and integrating over V we get 



B 4 + s(4 — s) 
—M s + (4 - s)aM s+1 + AaM s = ± >-M s 



(4 -a) (3 -a) s(s - I, , 
M s+2 -\ M, 



s-2 



s(4 - s) 



PiP, 



l-Ts-l 



o H o ■M)-«i— 2> 



(21) 



where Pq = Pxx/fnn, Pi = P xy /mn, and P2 = P yy /mn. In matrix notation, 



at 



-M. 



(22) 



where 



A ss >(a) 



4a(a) 



4 + s(4 - s) 



5 S / )S + (4 — s)a5 s ' !;j+ i 



(4- S )(3-s) 



s(s-l) 



s',s-2, 



(23) 



and is a bilinear combination of second-degree moments. The temporal evolution of the moments M s is governed by 
the eigenvalues X s (a) of the matrix A SyS >(a). At equilibrium (a = 0) the eigenvalues are Ao(0) = |, Ai(0) = A 2 (0) = 1, 
and A3(0) = A4(0) = |. At finite shear rate, the degeneracy is broken, (Ai, A 2 ) and (A3, A 4 ) becoming two conjugate 
pairs. The trace of the matrix A ss ' (a) is Ao(a)+2ReAi (a)+2ReA3(a) = ^+20a(a), which is a monotonically increasing 
function of the shear rate. On the other hand, as Fig. [j] shows, while the real parts of A1-A4 monotonically increase 
with the shear rate, the only real eigenvalue Ao decreases and eventually becomes negative for a > a c ~ 2.48553. This 
implies that for a > a c the fourth-degree moments diverge in the long-time limit. Of course, moments of a degree 
higher than four are also divergent for a > a c . As said in the Introduction, the existence of a critical sheattate a c 
was firstjobtained analytically in thepCase of Maxwell molecules with regular collision rates, both for d = 3ll3il3 and 
d = 2.E§EJ For instance, a c — 5.8470 in the case of two-dimensional Maxwell molecules with isotropic scattering. 
Thus, the singular behavior of the fourth-degree moments M s is not an artifact of the scattering model (|l6|). 

The divergence of M s for a > a c is a strong indication that the velocity distribution function reaches a steady 
state form that has an algebraic decay /(V) ~ I/~ 4 ~ CT ( a ) with a(a) < 2. Conversely, for a < a c one can expect that 
it (a) > 2, so that the fourth-degree moments converge, but moments of a higher degree 2k + 2 will diverge ..whenever 



the shear rate is such that a (a) < 2k. While this scenario has been supported by computer simulations,t2lll3 we are 
not aware of any derivation of this high-velocity tail from the Boltzmann equation, prior to the one presented in the 
next Section. 
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FIG. 1. Real parts of the eigenvalues X s associated with the fourth-degree moments, as functions of the shear rate. The 
arrow indicates the location of the critical shear rate a c ~ 2.48553 beyond which Ao < and, consequently, all the moments of 
degree four and higher diverge. 

IV. HIGH- VELOCITY TAIL 

From Eq. (|l^) we get the following equation for the steady state distribution function: 

-/• (24) 



8«V,A + ( i + 8a) ^ + £S 



P tj - 2p% d d 
ran dVi dVj 



This is still a complicated equation to solve for all velocities. However, we are interested here in the solution for large 
velocities. In that domain, the right-hand side of Eq. ( p4| ) is of order V~ 2 relative to the term C 2 f [cf. Eq. JlS|)], so 
it can be neglected. By the arguments given at the end of the previous Section, we look for self-consistent solutions 
with the asymptotic behavior 

f(V;a)~V- 4 -^<P(ip;a) (25) 

for large velocities. In Eq. ( p5| ) ip is the polar angle in velocity space, i.e., 

V x = V cos <yj, Vy = V sin ip, (26) 

$ is a function measuring the degree of anisotropy of the high-velocity distribution function, and a is the expo- 
nent characterizing the algebraic decay. These two quantities must be determined self-consistently. According to 



Eq. (|25D, the left-hand side of Eq. (24) is proportional to V~ 4 ~ a , while the right-hand side is proportional to V r ~ 6 ~ cr . 
Consequently, the latter must be neglected against the former in the limit of large velocities, as said before. 
In polar coordinates, the operators on the left-hand side of Eq. ( p^ ) become 

v 'W x =co " vsblvV w-" hl2 '<'^ (27) 



G 



(29) 



Thus, when Eq. 123) is inserted into Eq. (24) one gets the following linear second-order ordinary differential equation 



*"(¥>) - p [1 - cos(2<^)] $'(</>) - \p + 2gsin(2(^)] = 0, (30) 

where p = 4a, (3 = (2 + <r)(l + 8a), and q = 2a(4 + a). Equation ([30j) is a generalization of Mathieu's equatiorSlEI 
and reduces to it if p — 0. However, since in our case p ^ 0, we have to deal with Eq. (30) rather than with Mathieu's 
well-known equation. The symmetry property /(V) = /(—V) implies the periodicity condition <fr(ip) = $(92 + n). 
From a mathematical point of view, Eq. (|30| ) supports those periodic solutions provided that the parameter (3 takes a 
characteristic value (3(p, q). Since the shear rate a and the scaling coefficient a are related through Eqs. ( |l3| ) and (|l4|), 
only two of the three parameters (3, p, and q are independent. Thus, the characteristic value [3(p, q) translates into 
a(a). In case there are multiple solutions for a given a, the relevant solution is the one related to the dominant tail, 
which corresponds to the smallest value of a. In the subsequent mathematical analysis we consider /3, p, and q as 
independent parameters and only at the end we will take into account that in our physical problem (3 = (2 + a)(l + 8a), 
p = 4a, and q = 2a(4 + a). 

In order to find the periodic solutions to Eq. ( pp[ ) and the corresponding characteristic values, we write 

00 

$(<p)= J2 C ™z 2imv i C- m = C* m . (31) 

m— — 00 

If the above is substituted into Eq. ([50|), one obtains the recurrence relations 

C c m _! - 0, (32) 

where 

fJ- m = q-mp, ^ m = 2mp-i((3 + Am 2 ). (33) 

In particular, setting m = 0, 

C = -i^(C x -C7_ 1 )=2^ImC 1 . (34) 

Therefore, by applying Eq. ( |32] ) to m > 1, all the coefficients C m , m > 2, can be obtained recursively from Ci. The 
coefficients C m with m < —2 are then given by the symmetry relation C_ m = C* n . The characteristic values arise 
from the convergence condition 

km \D m \ < 1, D m = %±i. (35) 

The recurrence relation for the ratios D m is 

fi m+ iD m + 7 TO - /Lm+iD^j = 0. (36) 
Thus, all the coefficients D m with m > 1 are obtained from Do, the latter being subject to the compatibility condition 

ImDo ^ 2(h>y (37) 

as follows from Eq. ([34j). Equation ( |36|) shows that there are two possible asymptotic behaviors of D m for large m. 
Either 

4i 

D m « m (38) 

or 

D m « m" 1 . (39) 
4i 
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The latter is the only one consistent with the convergence condition (f55|). This allows us to use (36) to develop D n 
as the continued fraction 



D,, 



(J>-r 



7m+l + ^771+2^771+1 

M-m ^m+2/i-m-l Mm+3A*-Jn-2 



7m+l+ 7m+2+ 7777+3 + 

Finally, the characteristic values are obtained from the compatibility condition (|37|), namely 

Mo M3M-2 \ /? 



Im 



7i+ 72+ 73+ 



(40a) 
(40b) 

(41) 



Once the relationship /3(jp, g) is found, the ratios D m are obtained from Eq. (40b) or, equivalently, from Eq. ( |40q ) for 
m = and from Eq. (|36] ) for m > 1. The angular distribution function is then given by Eq. (|3l]) with 



777— 1 



m > 1, 



777' = 



where the particular value Co = 1/27T has been chosen to verify the normalization condition 

<-2ir 

dip$(ip) = 1. 



(42) 



(43) 



Since Eqs. ( |40b| ) and (|4l]) involve infinite continued fractions, in practice we proceed by setting = — {p/4i)N~ l 
for a certain large value of N and then using (40a) to get D m for m < N — 1 up to m = 0. We have typically taken 
N ~ 10 3 and have checked that the results are rather insensitive to a further increase of N. 

Before going back to our original physical problem and take into account the expressions of p, and q in terms of 
a and er, it is worthwhile exploring the following limiting cases. 

(i) Consider first that the ratio q/p is equal to an integer number, q/p = k + 2. This means that fik+2 — 0, so that 
Eq. (f4l|) becomes 



Im 



Mo M2M-1 M3M-2 Mfc+iM-fc 







,7i+ 72 + 73 + 7fc+i 7 2(fc + l)p' 

which gives rise to an algebraic equation of degree 2k + 3 for (3. For instance, 

P 3 + 8/3 2 + 16/3 - 16p 2 = 0, 



(44) 



(45) 



I3 b + 40/3 4 + 528/3 d + 16(160 - 2lp A )fi 2 + 256(16 - 15p A )(3 - 12288?/ = 0, 



(46) 



for k = and fc = 1, r espectively. Next, Eq. (40a) gives D m for m < k in a closed form, while the infinite 
continued fraction ( 40b ) must be used for m > k + 1. For instance, in the case k = 0, 



4p 



4p- (16 + ' 



4p - (16 + 0)i 



Yip 1 - (4 + 0) (16 + 0) - 6p(8 +0)i' 



(47) 



(ii) We now consider the limit q — > 00 with p = finite. In that case, as verified later by consistency, the characteristic 
value scales as (3 — > q(i\, so that \x m — ► g and 7 m — + —iqfii. Consequently, Eq. ( pq) becomes 



*A - d- 1 .! = 0. 



(48) 



The solution consistent with the convergence of the series and with Eq. (37) is simply D m = e 2iVo , where ipo 
is an angle defined by 



01 

sin 2(^0 = — ^ 



(49) 



8 



This implies that < j3\ < 2. The function $ is 

1 oo 



2w 

m —-co 



-[5(<p-ip ) + 5(<p-(p +TT)]. (50) 



(iii) As a third limiting case, let us consider q — > oo, p — > oo, g/p = finite. It is then expected that (3/q — > 0, so that 
7 m — » 2mp. The consistent solution to Eq. ( |36|) is Z? m = 1, which corresponds to 

*(¥>) = ~[%)+%-t)]. (51) 



Numerical analysis shows that in this limit the characteristic value (3 scales as (3 — (fcp 2 ) 1 ^ 3 , where 02 depends 
on the ratio q/p. In the special case of q/p = integer, the values of [3 2 can be obtained from the algebraic 
equation (Q). In particular, we get j3 2 = 16, f3 2 = 336, and f3 2 — 1008 ± 96\/79 for q/p = 2, q/p = 3, and 
q/p = 4, respectively. 

Now we apply the above analysis to the physical values (3 = (2 + cr)(l + 8a), p — 4a, and q = 2a(4 + a), where 
a and a are related by Eqs. (|l|) and (|i"i|). The case (i) above, namely q/p — k + 2 with k integer, corresponds to 
fj = 2k, so that the solution to Eq. (^4|) gives the critical value of the shear rate beyond which all the moments of a 
degree equal to or larger than 2k + 2 diverge. For k = 0, Eq. ( filil ) gives rise to a cubic equation in a with no positive 
real root. This is expected since, as seen in Sec. [n], the second-degree moments converge. For k — 1, Eq. (^6|) gives 
rise to the same quintic equation as obtained from det A ss / — 0, where the matrix A ss r is given by Eq. (p3[). The 
solution to this equation is a ~ 0.61796, which corresponds to a ~ 2.48553, i.e., the critical value a c found in Sec. || 
from the time evolution of the fourth-degree moments. In a similar way, we find the values a ~ 1.11175, a ~ 0.87611, 
a ~ 0.77052, and a ~ 0.70842 for a — 4, a — 6, a — 8, and a — 10, respectively. It is interesting to note that for a = 8 
and a = 1Q there exists a larger second solution (a ~ 52.5 and a ~ 7.3, respectively). This indicates that, in addition 
to the smallest eigenvalue, another of the eigenvalues governing the time evolution of the moments of degree 10 and 
12, respectively, becomes negative at those large shear rates. This property has been observed before in the case of 
more realistic scattering laws.Ej Of course, whenever there are different values of a consistent with a given value of 
cr, the actual function a(a) corresponds to the smallest value of a. Equivalently, if the characteristic equation ( |4l| ) 
gives more than a value of a for a given shear rate, the function <r(a) is defined by the smallest value. 

Now we consider the case (ii) above, which corresponds to the situation a — > oo. For which value of the shear 
rate does this happen? Mathematically speaking, there are an infinite number of solutions, namely the solutions to 
(1 + 8a)/2a = f3\ with < j3\ < 2. Taking into account Eq. (Jl3|), this is equivalent to the cubic equation 

8f3fa 3 + 4(5/3? - 16)a 2 + 6(3 x a -9 = 0. (52) 

The real solution to this equation is a monotonically decreasing function of (3\. Therefore, the physical value of the 
shear rate at which a — * oo corresponds to (3\ — 2. This defines the threshold shear rate ath = (v57 /96 + 67/864) 1 / 3 — 
(v / 57/96-67/864) 1 / 3 -l/12 ~ 0.352047 beyond which the high-velocity tail has the form (||). Conversely, if a < a th , 
the decay of the velocity distribution function is more rapid than any algebraic tail and all the velocity moments are 
convergent . ■ Thifi is in contrast with the results obtained from realistic scattering laws, which suggest that ath — > in 
those casesMIl 3 ] Since f3 x = 2, Eq. @ yields ip = 3tt/4 and Eq. @ becomes 

lim $(ip; a) = i 

Next, we investigate the limit of large shear rates, a — * oo. In that case, a w a 2 / 3 /2, on account of Eq. (|l3|). 
The question now is, what is the corresponding value er(oo) = cr m in? Computer simulations for Maxwell molecules 
interacting with an isotropic scattering law suggest that a m i n = 0. However, this is not the case with the scattering law 
(|l6|). According to the case (iii) above, the value of cr min must be such that (3 2 — 4(2 + er min ) 3 for q/p = 2 + cr min /2. 
A numerical calculation gives a ~ 1.252, i.e., (3 2 — 137.6 at q/p = 2.626, which is consistent with the fact that 
(3% = 16 and /?2 = 336 at q/p = 2 and q/p — 3, respectively. The minimum value <r m i n means that, no matter how 
large the shear rate is, not only the second-degree moments are finite, but so are moments of the form (V k ) with 
k < 2 + (Tmin — 3.252. On the other hand, the high-velocity distribution function becomes strongly anisotropic for 
large shear rates since, according to Eq. (|5l|), 



Ifi 7T 

^ 4 



lim <r(a) = oo. 



(53) 







lim $(p; a) = - [S(tp) + 8(<p - tt)] , lim a {a) = a min . (54) 

a — >oc Z & — >oo 



The complete dependence of a on the shear rate is shown in Fig. where the dotted lines indicate the locations 
of a t h and a m j n . The region near a = a t h is shown in Fig. || The results can be well fitted to the power law 
a w e~ 02 (a — a t h)~ 193 . The behavior for large shear rates is shown in Fig. [|, where the results can be fitted to 




FIG. 3. Log-log plot of a versus a — ath- The solid line is a linear fit. 
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FIG. 4. Log-log plot of a — a m i n versus a. The solid line is a linear fit. 
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In addition to the function a (a), the solution to the problem gives the angular distribution <j>(y>;a) for any shear 
rate a > ath- The results indicate that in the high-velocity domain the distribution is highly anisotropic, with the 
angular distribution sharply concentrated around a preferred orientation angle Ipioi) that rotates from tp = 37r/4, 7tp/4 
when a — ► a t h, Eq. (^3|), to tp = 0, tt when a — * a^, Eq. (p4]). This transition is illustrated in Fig. || The anisotropic 
behavior observed in Fig. || is consistent with the one shown in Fig. 8 of Ref. |l3| in the case of computer simulations 
of Maxwell molecules with a constant collision rate. 



330 




330 



FIG. 5. Polar diagram of $(y>) for (a) a = 0.354652 (a = 10 5 ), (b) a = 0.442229 {a = 10 2 ), (c) a = 0.708415 (a = 10), (d) 
a = 2.485530 (a = 2), (e) a = 10 {a = 1.427204), and (f) a = 10 2 (a = 1.282082). Case (a) corresponds to a shear rate slightly 
larger than the threshold value ath — 0.3520, while case (f) corresponds to an exponent a slightly larger than the minimum 
value (Tmin — 1.252. Note that the scales in cases (a), (e), and (f) are different from the scales in cases (b)-(d). 
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V. CONCLUDING REMARKS 



In this paper we have shown that the stationary solution of the (thermostatted) Boltzmann-Fokker-Planck equation 
for a two-dimensional gas of Maxwell molecules under uniform shear flow exhibits an algebraic high-velocity tail of 
the form given by Eq. (|2^), where a (a) is a decreasing function of the shear rate a. As a consequence, all the velocity 
moments of a degree equal to or larger than 2 + a (a) diverge. The word "thermostat" does not mean in the present 
context that an external drag force is necessarily applied on the system; its effect for Maxwell molecules is equivalent 
to a rescaling of the velocities with respect to the thermal velocity and this is the point of view adopted in this paper. 
The restriction to an idealized scattering model with grazing collisions, Eq. (|l^), allows us to replace the complicated 
structure of the Boltzmann collision operator by a much more manageable Fokker-Planck operator, Eq. (|T7|). This 
operator is formally linear in velocity space, but the original bilinear nature of the collision operator appears via 
the dependence of some of the coefficients on the distribution function through the pressure tensor. On the other- 
hand, those nonlinear terms are negligible in the limit of large velocities and thus one arrives at a linear second-order 
ordinary differential equation for Eq. (|30|). Despite the linearity of the problem posed by Eq. (p0|), it applies to 
states arbitrarily far from equilibrium. As a matter of fact, the nonlinear nature of the underlying state is still present 
through the dependence of the "thermostat" (or scaling) parameter a on the shear rate a, Eq. (jl4|). The structure of 
Eq. ( J30| ) is reminiscent of Mathieu's equation; in fact the latter is obtained by formally setting the parameter p = 0, 
although this is not possible in the physical problem at hand. The symmetry property /(V) = /(—V) translates 
into the periodicity condition $(<p) = $(<^ + n). Enforcements of this condition on the solutions to Eq. ( |3C| ) allows 
us to obtain a(a) from the solution of the continued fraction representation ([51), the latter becoming an algebraic 
equation when a = 2,4,6, .. .. The results show that cr{a) is lower bounded by the value er m i n ~ 1.252 (which is 
asymptotically reached in the limit a — > oo) and goes to infinity when the shear rate approaches the threshold value 
ath — 0.352 (in units of the collision frequency). The first property implies that all the moments of degree equal to or 
smaller than 2 + a m i n ~ 3.252 are finite in the steady state, regardless of how large the shear rate is. As a consequence 
of the second property, there exists a window of shear rates < a < ath where the high-velocity population decays 
more rapidly than algebraically (e.g., as a stretched exponential), all the moments being finite. The lobular shape of 
the orientation distribution $(y) [cf. Fig. |5|] shows that the high-velocity population is highly anisotropic, a feature 
already observed in Monte Carlo simulations for the isotropic scattering modelila Most of the particles having a large 
(peculiar) velocity V move along a narrow bunch of directions around a preferred direction characterized by the 
(equivalent) angles <J5(a) and <p(a) +7r. The angle <p(a) rotates counter-clockwise from tp — 37r/4 (i.e., V y /V x — —1) at 
a = ath to tp = 7r (i.e., V y /V x — > _ ) in the limit a — > oo. The width of the angular distribution is zero in both limits, 
being maximum at a shear rate awl. 

It is worth remarking that, despite the simplicity and artificiality of the scattering model considered in this paper, 
it succeeds in capturing the most relevant features that were expected-Xin the basis of Monte Carlo simulations and 
moment method results in the case of more realistic scattering laws:E3E3 the existence of a high-velocity tail with a 
monotonically decreasing exponent <r(a) and with a strongly anisotropic orientation distribution $(<y9;a). Thus, the 
asymptotic behavior (^5|) is not an artifact of the specific model (|l7|), but a general property of Maxwell molecules 
under uniform shear flow. Some other features, however, are possibly peculiar of the singular scattering law (|l6|). For 
instance, the results of Ref. [l3| for the isotropic scattering model seem to indicate that a(a) only diverges when a — > 
(i.e., ath = 0) and that lim a ^ 00 CT(a) = (i.e., cr m j n = 0). It is possible that the extreme anisotropy of the scattering 
model ( |l6| ) makes the high-velocity tail phenomenon to be milder than in the general case, thus yielding non-zero 
values for ath and er m i n . On a different vein, it must be said that since the problem reduces to a linear equation, we 
have not determined either the coefficient measuring the amplitude of the tail (^5|) or ±he order of magnitude of the 
characteristic velocity c, such that (|25| ) applies when V ^> c. Monte Carlo simulations^ show that the amplitude is 
a decreasing function of the shear rate, while c is rather insensitive to the shear rate. The value of c seems to be 
of the order of the thermal velocity vq = ^JlksT jm, so that the asymptotic behavior (^B|) is reached in practice for 
V > IOuq. The theoretical confirmation of both properties would require the analysis of the full kinetic equation (p3) 
and this is beyond the scope of this work. 

It must be remarked that in this paper we have considered the replacement ( |l7| ) as a collision model by itself. On 
the other hand, it is known that the right-hand side of ([l7]) actually represents the contribution to the true collision 
operator associated with grazing collisions.El Therefore, it is tempting to conjecture that, in general, the high-velocity 
tail is produced. by grazing collisions of high- velocity particles. Since the expansion leading to (jnj) is non-uniform for 
high velocities,El the analysis of this paper does not provide a rigorous proof of the above conjecture but, at most, a 
strong indication of it. 

The results presented in this paper, along with those of Refs. |Io|-|l4|, show that the population of high energy 
levels in the uniform shear flow, and possibly in many other nonequilibrium states, is far greater than in the corre- 
sponding equilibrium state described by the Maxwell-Boltzmann distribution function. Among other examples of 
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nonequilibcium states exhibiting algebraic high- velocity tails we can mention the-viscous longitudinal flow of Maxwell 
moleculesc3 and the homogeneous cooling state of inelastic Maxwell molecules.E3 As a consequence of these nonequi- 
librium overpopulation effects, many more particles than at equilibrium can be available for surmounting the energy 
barriers of chemical and nuclear reactions which, consequently, would proceed faster. We also guess that fluctuations 
are enhanced with respect to the equilibrium ones, but the analysis of this expectation would require the treatment 
of the fluctuating Boltzmann equation.EH This is an interesting avenue for further theoretical work. 
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APPENDIX A: DERIVATION OF EQ. (|ll|) 

The Boltzmann collision operator for Maxwell molecules is 

Q[Vi|/,/] = J dv 2 J aiK(0)[f(yi)f(v' 2 )-f(Vi)f(y2)], (Ai) 

where 

Vi=Vi-(g-&)$, V 2 =V 2 + (g-£)a. (A2) 

Here, g = Vi — V 2 is the relative velocity and <r is a unit vector lying on the scattering plane and making an angle 
equal to 5(71" — 6) with the vector g. 

Given an arbitrary function if (Vi), standard manipulations yield 

J dViff(V 1 )Q[V 1 |/,/] = l J dVi J dV 2 J dnK(e)f(V 1 )f(V 2 )[H(V[) + H(V' 2 )-H(V 1 )-H(V 2 )}. (A3) 

In the particular case of if(Vi) = mViVx, 

H(V[) + H(V' 2 ) - H (Vi) - if (V 2 ) = m (g • a) [2 (g ffg - gS] . (A4) 

Taking into account the identities 

J dnK{6)(g-*)a = B lg , (A5) 

J dnK(9) (g • af aa = i? 2 gg - ^ ~ f 2 (gg - .9 2 l) , (A6) 
where I is the d x d unit tensor and 

B k = J dQK(e) sm 2k ^, (A7) 

we finally have 

J dVt mVxViQlVil/,/] - -jL^Bi-Bi) J dVi J dV 2 f(V 1 )f(V 2 )m (gg- ~g 2 \) 

- d (B 1 -B 2 )n(P-p\). (A8) 



d-1 

This allows us to identify the effective collision frequency of Eq. (0) as v = n{B\ — B 2 )2d/ (d — 1) 
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APPENDIX B: DERIVATION OF EQ. (17) 



Let us start rewriting Eq. (A3) as 

I[H] = J dV x H{Y 1 )Q[Y 1 \f,f] = JdVx J dV 2 J dQK(6)f(V 1 )f(V 2 ) [H(V[) - H(V X )] , (Bl) 



where H(Vi) is an arbitrary function. Now, according to the scattering law (|16[), (g • a) — gsin(e/2), so that 
— Vi ~ e. This justifies the approximation 



ff(Vi) - H(V X ) 



aff(Vi) id 2 H(v 1 ) , , 

— {g-cr)cri + - 



dV u 



2 dVudV-i 



(g • a) didj 



Consequently, 



dH(Vi) ld 2 H(V 1 ) 



-9i + 7: 



2 0Vii0V y 



n Bl — B 2 , Of \1 1 



(B2) 



(B3) 



where we have made use of (A.5) and flAq ). On the other hand, lim e ^o B 2 = 0, while B\ — > (v/n)(d — l)/2d remains 
constant. Integrating over V2 in Eq. (|B3|), we have 



I[H] 



d-1 
2d 



v \ dV 1 H(V 1 ) 



d 



dV u 



d 



2(d - 1) dVj 



ran 



/(Vi), 



where we have integrated by parts. Since i/(Vi) is arbitrary, it follows that 



Q[Vi|/,J] 



d-l d 
-v- 



4d dV u 



1 



d — 1 \ mn 



/(Vi), 



(B4) 



(B5) 



1 J. W. Dufty, J. J. Brey, and A. Santos, in Molecular-Dynamics Simulation of Statistical-Mechanical Systems, edited by G. 
Ciccotti and W. G. Hoover (North-Holland, Amsterdam, 1986), pp. 294-303. 

2 A. W. Lees and S. F. Edwards, J. Phys. C 5, 1921 (1972). 

3 J. M. Montanero, V. Garzo, and A. Santos, Phys. Fluids 12, 3060 (2001). 

4 A. Santos, V. Garzo, and J. J. Brey, Phys. Rev. A 46, 8018 (1992). 

J C. Cercignani, Theory and Application of the Boltzmann Equation (Scottish Academic Press, Edinburgh, 1975). 

6 C. Cercignani, Mathematical Methods in Kinetic Theory (Plenum Press, New York, 1990). 

7 E. Ikenberry and C. Truesdell, J. Rat. Mech. Anal. 5, 1 (1956). 

8 C. Truesdell, J. Rat. Mech. Anal. 5, 55 (1956). 

9 C. Truesdell and R. G. Muncaster, Fundamentals of Maxwell's Kinetic Theory of a Simple Monatomic Gas (Academic Press, 
New York, 1980). 

10 A. Santos and V. Garzo, Physica A 213, 409 (1995). 

11 A. Santos, V. Garzo, J. J. Brey, and J. W. Dufty, Phys. Rev. Lett. 71, 3971 (1993); 72, 1392 (Erratum) (1994). 

12 J. M. Montanero, A. Santos, and V. Garzo, Phys. Rev. E 53, 1269 (1996). 

13 J. M. Montanero, A. Santos, and V. Garzo, J. Stat. Phys. 88, 1165 (1997). 

14 J. M. Montanero, V. Garzo, and A. Santos, in Rarefied Gas Dynamics 20 th , edited by C. Shen (Peking University Press, 
Beijing, 1997), pp. 113-117. 

15 J. R. Dorfman and H. van Beijeren, in Statistical Mechanics, Part B, edited by B. J. Berne (Plenum, New York, 1977), pp. 
65-179. 

16 J. W. Dufty, A. Santos, J. J. Brey, and R. F. Rodriguez, Phys. Rev. A 33, 459 (1986). 

17 C. Cercignani, Transport Theory Stat. Phys. 182, 125 (1989). 

18 C. Cercignani and S. Cortesse, J. Stat. Phys. 75, 817 (1994). 

19 M. H. Ernst, Phys. Rep. 78, 1 (1981). 

20 A. V. Bobylev, Sov. Sci. Rev. C Math. Phys. 7, 111 (1988). 

21 J. A. McLennan, Introduction to N on- equilibrium Statistical Mechanics (Prentice Hall, Englewood Cliffs, N. J., 1989). 



14 



22 J. M. Montanero, A. Santos, and V. Garzo, Phys. Fluids 8, 1981 (1996). 

23 W. G. Hoover, Ann. Rev. Phys. Chem. 34, 103 (1983). 

24 D. J. Evans and G. P. Morriss, Comput. Phys. Rep. 1, 299 (1984). 

25 D. J. Evans and G. P. Morriss, Statistical Mechanics of Nonequilibrium Liquids (Academic Press, London, 1990). 

26 Handbook of Mathematical Functions, edited by M. Abramowitz and I. A. Stegun (Dover, New York, 1972), Ch. 20. 

Table of Integrals, Series and Products, edited by I. S. Gradshteyn and I. M. Ryzhik (Academic Press, New York, 1965), Ch. 
8.6. 

28 A. Santos, Phys. Rev. E 62, 6597 (2000) and |cond-mat /0002240 . 

29 P. L. Krapivsky and E. Ben-Nairn, J. Phys. A: Math. Gen. 35, L147 (2002); A. Baldassarri, U. Marini Bettolo Marconi, and 
A. Puglisi, Europhys. Lett. 58, 14 (2002); M. H. Ernst and R. Brito, Europhys. Lett. 58, 182 (2002). 

30 E. G. D. Cohen, in Fundamental Problems in Statistical Mechanics V, edited by E. G. D. Cohen (North-Holland, Amsterdam, 
1980), pp. 235-248. 



15 



